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Abstract 

We use gauge-gravity duality to study the temperature dependence of the zero sound mode and 
the fundamental matter diffusion mode in the strongly coupled = 4 SU{Nc) supersymmetric 
Yang- Mills theory with Nf J\f = 2 hypermultiplets in the Nc ^ I, Nc ^ Nf limit, which is 
holographically realized via the D3/D7 brane system. In the high density limit /i ^ T, three 
regimes can be identified in the behavior of these modes, analogous to the collisionless quantum, 
collisionless thermal and hydrodynamic regimes of a Landau Fermi-liquid. The transitions between 
the three regimes are characterized by the parameters r//i and (T/^u)^ respectively, and in each 
of these regimes the modes have a distinctively different temperature and momentum dependence. 
The collisionless-hydrodynamic transition occurs when the zero sound poles of the density-density 
correlator in the complex frequency plane collide on the imaginary axis to produce a hydrodynamic 
diffusion pole. We observe that the properties characteristic of a Landau Fermi-liquid zero sound 
mode are present in the D3/D7 system despite the atypical T^//i'^ temperature scaling of the 
specific heat and an apparent lack of a directly identifiable Fermi surface. 
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I. INTRODUCTION 



The AdS/CFT correspondence [Mi^ ^^d, more generally, gauge-gravity duality have been 
used extensively for studies of the thermodynamics and transport properties of strongly- 
interacting quantum field theories at finite temperature and density. Applications of the 
correspondence to real-world physics initially concentrated on QCD and, in particular, the 
quark-gluon plasma s], 6|, but in recent years there has been a surge of interest in using these 
tools to study the physics of condensed matter systems (see e.g. {t-Q] for an introduction 
to the field). 

In this paper we investigate the behaviour of the collective modes of the = 4 SU (Nc) 
supersymmetric Yang-Mills theory coupled to Nf M = 2 fundamental hypermultiplets at 
infinitely large 't Hooft coupling A = Qym^c in the limit ^ 1., ^ Nf and at finite 
temperature T and finite density of the fundamental matter d. The holographic dual of this 
theory is provided by embedding Nf D7-branes in the gravitational background created by 
Nc D3-branes and treating the D7-branes as probes \u^. Although a seemingly elaborate 
construction, this is one of the simplest known finite-density field theories for which an 
explicit dual is known, and it has some interesting properties. 

At strictly zero temperature and zero hypermultiplet mass, the theory supports a col- 
lective excitation which appears as a pole in the density-density correlator Such an 
excitation, found in [Uj using the dual gravity methods and thus referred to here as the 
"holographic zero sound" , is reminiscent of the zero sound mode predicted to exist by Landau 



in a class of 
Helium-3 



18 



ermi liquid systems 



12| (see also |l3l-ll7]) and subsequently observed in liquid 



19| . In Landau Fermi-liquids (LFLs), the zero sound mode arises due to (non- 



thermal) interactions between the constituent fermions which result in oscillations of the 
Fermi surface. The holographic zero sound mode in the D3/D7 system at zero temperature 
has a speed equal to that of the ordinary (first) sound and an attentuation proportional to 
the square of the momentum. This is identical to what one finds in the Fermi liquid models 
where the interaction strength (parameterised by the Fermi liquid coefficient Fq) approaches 
infinity [16]. At the same time, the heat capacity of the D3/D7 system is proportional to 
at low temperatures whereas in normal Fermi liquids it is proportional to T. 

Similar investigations have been made in other string-theoretic constructions - in 20|] it 
was shown that the holographic zero sound mode persists at T = when the hypermultiplet 
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FIG. 1. Relative scales in the hydrodynamic, collisionless thermal and collisionless quantum regimes 
of a Landau Fermi-liquid. 



is given a finite mass, in 2l|] a similar mode was found in the T = D4/D8/D8 theory at 



finite density, in 



22| a "zero sound" mode was reported to exist in the (l+l)-dimensional 



theory on the D3/D3 intersection and in 23| a low temperature zero sound mode was found 
in the theory on the (2+l)-dimensional D3/D7 brane intersection. In the D4/D8/D8 case - 
where the only fundamental matter present is fermions - the heat capacity is proportional 
to T as expected for a Landau Fermi-liquid. However, the imaginary part of the zero sound 
mode in this case has an unconventional dependence, which is at odds with the predictions 
of LFL theory. Holographic zero sound modes have also been discovered in theories con- 



sisting of probe DBI actions embedded in Lifshitz (i.e. non-relativistic) spacetimes [24l. |25|. 
Such modes have a density-dependent speed, and their properties depend upon the critical 
exponent of the spacetime. In addition to these, the field theory dual to (3+l)-dimensional 
Einstein-Maxwell gravity with a cosmological constant was found to support a long-lived 
sound mode at zero temperature [26 1. 

A natural question to ask is what happens to the holographic zero sound mode when the 
temperature is turned on. In 27(], it was argued that at any infinitesimal temperature the 
zero sound mode of the D3/D7 theory is no longer dominant and instead density transport 
occurs mainly via diffusion. In this paper we show that this is not the case - at low tempera- 
tures, the holographic zero sound mode still dominates the density-density correlator with a 
dispersion relation very similar to that at T = 0. Furthermore, we show (numerically) that 
the real part of this mode is independent of temperature whereas its imaginary part receives 
corrections proportional to T^. At a sufficiently high temperature, there is a crossover to 



a hydro dynamic regime where the diffusion mode of 27|] dominates (such crossovers are 



common - see 28l433l| for other holographic examples of this). This behaviour is familiar as 
it is the same as that predicted by Landau's theory of Fermi liquids, to which we now turn 
to help clarify these results and to identify the relevant scales involved. 

TABLE I. Relative scales in a Landau Fermi-liquid 
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In Landau's theory of Fermi liquids, one assumes that the ground state is a degener- 
ate system of interacting fermionic quasiparticles. The theory describes small fluctuations 
around this ground state due to quasiparticle-quasihole excitations and/or collective excita- 
tions, and remains valid as long as these fluctuations are sufficiently small. The applicability 
conditions of Landau Fermi-liquid theory are given by the inequalities 



(1) 



where /i ~ d^f^ is the chemical potential. The first inequality guarantees that we are 
considering a degenerate (i.e. quantum) liquid, while the second ensures that the excitations 
are sufficiently macroscopic, i.e. their wavelengths are much larger than the characteristic 
interparticle distance (this condition also implies that the quasiparticles remain sufficiently 
close to the Fermi-surface: u ~ \ep — n\ <^ ji). The results of the Landau theory can be 
viewed as the leading order term in an expansion in powers of u/fi Q, 16]. 

In Fermi liquids, the zero sound mode is a longitudinal, gapless, collective excitation 
corresponding to oscillations of the Fermi surface around its equilibrium shape at zero tem- 
perature. Its dispersion relation u{q) = Vgq — iT^ contains a non-zero damping due to the 
mode's decay into quasiparticle-quasihole pairs. It is often more convenient to work with real 
frequency and complex wave vector: q{u) = u/vg + iTg. The damping rate is conveniently 
characterised by considering the argument of q, aigq{u)) = Img/Reg, as a function of the 
frequency. 

The zero sound mode persists at small, non-zero temperatures as well, but its properties 
are altered by the thermal collisions of the quasiparticles. Such collisions have a character- 



arg(q/Ai) 



A 




(T///)2 



FIG. 2. A sketch of the dependence of the sound mode damping on frequency in the hydrodynamic 
(I), colhsionless thermal (II) and collisionless quantum (III) regimes of a Landau Fermi-liquid. 
First sound propagates in region I while the zero sound mode exists in regions II and III. 

istic frequency u ~ 1/r, where r is the mean time between quasiparticle collisions. One 
distinguishes three regimes: the hydrodynamic regime, characterized by a; ^ z/, the col- 
lisionless thermal (classical) regime, with u ^ u, u T, and the collisionless quantum 
regime, w ^ z/, ^ T. The zero temperature zero sound mode persists essentially un- 
altered in the collisionless quantum regime, where thermal excitations are too weak and 
infrequent to influence it. As the temperature is increased, however, thermal excitations 
change the attenuation of the zero sound mode giving it nontrivial temperature dependence. 
In the hydrodynamic regime, thermal excitations destroy the zero sound mode completely. 
However, these excitations support the ordinary hydrodynamic sound mode with viscous 
damping, and the thermal diffusion mode. 

In the collisionless regime, the frequency u can be computed from kinetic theory applied 



where ep ^ fi. The decay rate of the zero sound mode in the collisionless regime is Fg ~ 
[igI . In the quantum collisionless regime one has then F^ ~ i^^//^ and argg(w) ^ lo/h 
whereas in the thermal collisionless regime Tq ~ T^//i and argg((X') ~ {T / ^)'^ fi / u . Finally, 
in the hydrodynamic regime, the conditions u <ti u and oj <^ T lead to u/fi (T///)^ 



to Fermi liquids [13|, [14 , 




(2) 



1/ ~ 
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(the "quantum" limit w 3> T, w <^ z/ lies outside of the LFL applicability range ([T])). The 
attenuation of hydrodynamic (first) sound is determined by the viscosity and is proportional 

to wVt^ yiiiliel- 



The dimensionless variables which are most convenient for identifying the three regimes 
are w//! and g//i. In these variables, the regions I, II and III corresponding to the hydrody- 
namic, coUisionless thermal and collisionless quantum regimes respectively, are separated by 
the scales (T/ /i)^ and T / fi. Alternatively, in the language of the "traditional" hydrodynamic 
variables u/T and q/T, the relevant scales are T//i, 1 and fi/T (see Table [land Fig. [1]). 

TABLE II. Sound attenuation coefficients in a Landau Fermi-liquid 





To, 




Argq 


Hydrodynamic regime 








M 


Collisionless thermal regime 








E 

LU 


Collisionless quantum regime 






M 





The sound attenuation constants in various regimes are shown in Table [TTl and Fig. |2j 
The temperature dependence of the sound attenuation coefficient F, is shown in Fig. [3l We 
shall use this information cLS cL SU. ggestive guide in our investigation of the holographic zero 
sound at finite temperature. 

In the massless D3/D7 system at low temperature, the chemical potential is proportional 
to the cubic root of the volume density d of the U{1) "baryon" chargcl Hi 

T 



a t 



(3) 



where a can be expressed using the Euler beta-function, a = 5(1/3, 7/6)/2 ^ 1.402. We 
shall study the D3/D7 theory in the limit T -C d^^^, to <^ d^^^ formally corresponding to 
the applicabihty regime ([T]) of Landau Fermi-liquid theory. The appropriate dimensionless 
variables are 

- q 7 d , , 

w = ^, q = -;T7^, d = j-^ . (4) 



A priori, we do not expect to find agreement with the LFL results outlined above since the 
D3/D7 system appears to be microscopically rather different, with no obviously detectable 



^ An explicit expression for the charge density operator involving fundamental fermions and complex scalars 
of the Af = 2 hypermultiplet is given in Appendix A of 3J| . 



6 
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FIG. 3. The temperature dependence of the sound attenuation coefficient Tg in various regimes of 
a Landau Fermi-liquid. Above: A sketch of the dependence in the hydrodynamic (I), colhsionless 
thermal (II) and colhsionless quantum (III) regimes. Below: Temperature dependence of the 
acoustic attenuation in liquid ^He at P=32 kPa measured at both 15.4 MHz (Q) and 45.5 MHz (□). 
The lines through the data correspond to logFg ~ 21ogT and logF^ ~ — 21ogT in the colhsionless 
thermal and hydrodynamic regimes, respectively, in agreement with Table [Hi Reprinted with 
permission from Abel et al. [18]. Copyright (1966) by the American Physical Society. 
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Fermi surface or long-lived quasiparticles in its vicinity. Nevertheless, we do find that the 
behaviour of the zero sound mode at finite temperature is qualitatively similar to the one 
predicted (and observed) in a Landau Fermi-liquid: 

i) The three regimes (hydrodynamic, collisionless thermal and collisionless quantum) are 
readily identified by analysing the behaviour of the lowest quasinormal frequencies and the 
spectral function of the charge density correlator. The hydrodynamic - collisionless thermal 
transition occurs at cue? ^/'^ ~ 1 and is most spectacularly manifested in the motion of the 
zero sound poles in the complex frequency plane: as the temperature is increased, the two 
poles (corresponding to zero sound propagation with velocities +Vs and —Vg in the field 
theory) recede into the complex plane, approximately tracing a circle, until they collide on 
the imaginary axis and form two poles with zero real part - one of these new poles recedes 
even deeper into the complex plane while the other approaches the origin. The latter pol e 
is the charge density diffusion mode (with the diffusion constant computed in [27 



characteristic of the hydrodynamic regime. Such behaviour was previously observed in 



31| 



and also in a (2+l)-dimensional holographic field theory in [23|], where it was correctly 
identified as the hydrodynamic - collisionless transition involving the zero sound mode. A 
similar transition between propagating and diffusive modes has been seen in a model of a 



holographic superconductor [33||. The second transition, between the collisionless thermal 
and the collisionless quantum regimes, is observed at 6j d^^^ ~ 1. 

As our investigation is limited to the current-current correlators, we are not able to follow 
the emergence of first sound in the hydrodynamic regime - this should appear as a pole in 
the energy-momentum tensor correlators of the field theory which are decoupled from the 
current-current correlators in the probe brane limit. Accordingly, we do not expect the 
details of the hydrodynamic - collisionless thermal transition described above to survive 
beyond the probe brane approximation: it seems more likely that the acoustic poles, rather 
than colliding on the imaginary axis as the temperature is raised, will recede deeper into 
the complex plane, reflecting the behavior of the attenuation curve in Fig. [3l and then come 
back close to the real axis again as the hydrodynamic sound poles. 

ii) The D3/D7 zero sound attenuation coefficients exhibit the momentum, temperature 
and density dependence typical of a Landau Fermi- liquid as shown in Table [III We find that 
the dependence of the acoustic damping upon frequency and temperature is qualitatively 
the same as shown in Figs. [2|, [3]in regions II and III. 
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iii) These results remain valid in the case of a non-vanishing hypermultiplet masso 
The structure of the paper is as follows. In Section [TTl we give a brief description of 
the gravitational dual of the D3/D7 field theory and the relevant, known properties of its 
Green's functions and spectral functions. In Section [XTTl we present our numerical results for 
the density-density spectral function (and the dominant pole of the corresponding Green's 
function) when T <^ d^. We identify the three regimes similar to those of a Landau Fermi- 
liquid, and describe in detail the behaviour of the collective modes as the temperature of 
the system is varied. We summarize our results and discuss how they may generalize to 
other holographic finite density systems in Section IIVI Some details relevant for the case of 
a massive hypermultiplet are relegated to appendices: Appendix |A] contains the action and 
equations of motion for the fluctuations, and Appendix [B] provides a derivation of the zero 
sound attenuation constant at zero temperature. 



II. THE D3/D7 SYSTEM 

The specific field theory whose elementary excitations at finite temperature and density 
we wish to investigate is (3+l)-dimensional Af = 4 SU{Nc) supersymmetric Yang-Mills 
theory coupled to Nf M = 2 fundamental hypermultiplets with a global U{Nf) flavour 
symmetry. It arises as the low-energy theory on the worldvolume of a set of Nc D3-branes and 
Nf D7-branes intersecting along (3-1-1 )-dimensions. Taking Nc oo with both A = Qym^c 
and Nf/Nc fixed, and subsequently taking A — )■ oo and Nf/Nc — ?■ 0, we obtain a classical 



gravitational dual to this field theory [10|, l36| : 

S 'S'adjoint ~l~ •^'fundamental ; 



•S'adjoint - NfToj j rf^^A/-det {gab + Fah), 



(5) 



where ^adjoint is the ten-dimensional supergravity action and T^j is the tension of a D7- 
brane. In this probe brane limit, the metric is fixed and it is 5'fundamentai which contains 
the dynamical information. For a zero temperature field theory, the contribution of the 
fundamental matter is the DBI action of Nf probe D7-branes extended along an AdS^ x 5*^ 
section of the (fixed) AdS^ x background spacetime generated by the D3-branes. In 



Within the region of thermodynamic stabihty of the theory. 
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Eq. ([5]), gab denotes the induced worldvolume metric on the D7-brane and Fab is the field 
strength of a worldvolume U{1) C U{Nf) gauge field. 

When the field theory is at a non-zero temperature, the background spacetime is that of 
an AdS-Schwarzschild black brane (with a horizon at r = r^) times a five-sphere: 



ds' - — 



l_UL]dt^ + dx' 



l-'-^V'^dr^ 



R^ds% 



(6) 



The D7-brane wraps an asymptotically AdS^ x section of the metric, with the horizon 
radius of the background related to the temperature of the field theory via T = r/z/vri?^. 
In terms of the dimensionless radial coordinate u = rfj/r"^, the metric can be written 

as 

^^fdt^ + dx^) + -^du" + {de^ + sin^ eds% + cos^ ^^43) , (7) 
^. In these coordinates, the horizon is located at m = 1 and the 



dn^ 



U 



where f{u) = \ 

boundary at m = 0. In equilibrium, the D7-brane embedding can be characterised by a 
single embedding coordinate 0[u) (which determines which S"^ section of the background 5*^ 
it wraps). The gauge field on the brane is dual to a global fiavor current in the field theory 
and thus turning on the time component of a f/(l) C U{Nf) gauge field At{u) on the brane 
corresponds to introducing a finite density d of the U{1) "baryon" charge in the field theory. 
In this case, it corresponds to a net density of fundamental fermions and scalars. 

The equations of motion for the background fields are obtained from the DBI action 

cos^ 6 



Si 



fundamental 



Nrj, 



dud X- 



1 + Au^fO'^ - A^A 



12 



(8) 



H 



where N = NfTmVsa is a normalization constant determined by the gauge-gravity duality 
dictionary 3^ and primes denote derivatives with respect to u. One of the equations of 



motion derived from the action ([8]) reduces to 




cos^ 6{u) + d'^u^ 



(9) 



where d = dR^/r'fj = d/ (ttT)^ is a dimensionless parameter of the field theory related to 



the net number density of "quarks" Ug in the field theory 



via 



d 



25/2 



y/XNfNcT^ ' 



(10) 



^ See footnote [TJ Note that our normalization of d is different from the one used in 
whereas in our case d ^ Uq/^/XNcN f . 



341: in 



34|, d 
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For a given density d, there is a corresponding dimensionless chemical potential given by 

(11) 



d 



^=2 



du\ 



'i + 4/(m)m2^'(m) 



cos^ 9{u) + d^u^ 



which is related to the field theory chemical potential /li^T by /i = uhiT = y f In the 
massless case, the integral in Eq. (ITT]) can be computed exactly [ll|: it reduces to (E]) in 



the low temperature limit. The important scaling (not entirely o 



generically accompanied by mass and temperature corrections 
The equation of motion for the embedding coordinate is 



3vious from Eq. ( ITT]) and 
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37( 1) to emphasize is 



T 




f{u) cos^ 9{u)e'{u) \ 3 cos^ e{u) sin e{u) ^/G{ 



u 




u) 



where 



G{u) = cos^ e{u) 



4^3 



(12) 



(13) 



cos^ 6(u) + u^d"^ _ 
Eq. ( 1T3|) has no known generic analytic solution and must be solved numericalljo. Near the 

boundary, the solution has the form 

rh 



e{u) 



(14) 



where rh is related to the bare masqj of the fundamental matter in units of temperature via 



39| 



2 M„ 

^ = ^^- (15) 

To study the temperature dependence of the theory at finite density, it will be more conve- 
nient to use the field theory mass normalized with respect to density rather than tempera- 
ture: 



m 



m 
d^ 



2tt Ma 



\/A di 

Note that 9{u) = for a massless hypermultiplet, and thus the dual description greatly 
simplifies in this case - there is only one non-trivial background field At{u) which has the 
simple equation of motion 



rnd 



d\^ 



-1/2 



(17) 



^ A method to approximate the sohition in the low temperature Hmit has been proposed in 
^ The exphcit field theory "quark" mass operator is given in Appendix A of [34i]. 
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The thermodynamics of the system can be determined by computing the on-shell ac- 
tion. In the massless case, the specific heat schematically has the following temperature 
dependence 

cv ~ NlT"" + • • ■ + XNfN.ryd +■■■ , (18) 

where the first term comes from the adjoint A/" = 4 SYM degrees of freedom. The pressure 
P, the entropy density s and the ener gy d ensity e can be computed from the thermodynamic 
potential f2(/i,T) given explicitly in llj: the entropy density is finite in the T — )■ limit, 



s ~ yU^, and the equation of state is e = 3P for all T and fi (and thus the speed of first sound 
is Vs = l/y/S). In certain regimes, the system becomes thermodynamically unstable: there 



appears to be an instability at 
and large mass m ^ 1 , d ^ 1 



ow density d < 



0.00315 



32 



34| . and at high density 



In addition, when d ^ and m ^ there is a lower 



limit on the possible value of u, below which there is a phase transition to a Minkowski 



embedding with d = 



21, 



40l443l|. The values of the parameters considered in this paper 



are outside of these regimes. 

To determine the field theory excitations, one considers fiuctuations of the background 
bulk fields 

6{u) 6{u) + (j){u,z,t) , 

At{u) ^ At{u) + at{u,z,t), (19) 
Ai{u) ai{u, z,t) , 

i = x,y, z, where we have chosen the fiuctuations to depend on time, the radial coordinate 
and one of the spatial coordinates (z) only (the latter is possible due to the isotropy of the 
theory). Introducing Fourier components for the embedding fiuctuation 

du dq 



{u,z,t) 



(f){u,uj,q). 



(20) 



and similarly for a^{u,z,t), and expanding the DBI action (jH]) to quadratic order in fiuc- 
tuations, one finds the resulting action and the corresponding equations of motion. The 
transverse modes dx, dy decouple from the longitudinal modes a^, dz and 0. In the follow- 
ing, we shall focus on the longitudinal modes as they are the ones that encode the fate of 
zero sound at finite temperature. 

For a non-zero hypermultiplet mass, the equations of motion lead to a pair of coupled 
differential equations for the embedding fluctuation (j){u,u,q) and the gauge- invariant com- 
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bination (c.f. 44| ) 

Z{u,u,q) = ud^{u,u,q) + qat{u,u,q) , (21) 

where the dimensionless variables u and q were introduced in Eq. The equations of 
motion and the action for the longitudinal fluctuations are given in Appendix |A1 

In the zero mass limit, the equations for these two modes decouple and the gauge-invariant 
combination (12T]) obeys the equation of motion 



d 

du 



f{u) Z' 



g[u) (u;2 



+ 



daZ 



0, 



(22) 



where g{u) = (1 + d u ) . In this limit, the longitudinal part of the on-shell action is 



c(2) 
'-'long. 



dujdq f Z{u, —u, —q)Z'{u, u, q) 



(27r) 



^ (a;2 - q^fg) 



Ub 



(23) 



In this massless case, we obtain the longitudinal retarded Green's functions from the 
gravitational fields via the usual procedure 45 1 

oj"^ Z'{e,U!,q) 



G%j. (u;,g) = -lim2iVr^ 



f Z{e,u,q) ' 



(24) 



where the — o)^ factor comes from the definition (I2T]) . and Z{u,uj,q) is the solution obeying 
ingoing boundary conditions at the horizon: Z{u,uj,q) ~ (1 — uY as m — )• 1, where 7 = 
—%{jjd^l'^ j 4. Poles of the retarded Green's function are determined by the values of u) (q) for 
which the solution obeying the ingoing condition at the horizon vanishes at the boundary 



44 



451 ] ■ The density-density correlation function follows trivially from the conservation 



of current Gjtjt {00, q) = q"^ Gj^jz/u"^ and has the same poles as the longitudinal Green's 
function ( 12^ . The spectral functions of these operators are then obtained by taking the 
imaginary part: 



Xzz (w, q) = -2 Im [G%j. {u, q)] , Xtt (w, q) = -2 Im [G^tjt {u, q 



(25) 



with Xtt = Xzz/^'^- At zero temperature and density, the form of the longitudinal spectral 
function is known analytically 46| : 

NfN, 



Xzz (w, q) 



4:71 



{u^ - g2) e {u' - q') sgn 



(26) 



where is the Heaviside step function. 



13 



In the m 7^ case, the coupled equations of motion for bulk fluctuations imply that the 
dual field theory operators mix and the method to determine the retarded Green's functions 
is more involved (see Appendix [X] for the details). 

At zero mass, high densities co,q <ti and strictly zero ternperature, the dominant pole 
of the correlators Gjtjt and Gjz jz has the dispersion relation [ill, 



- I ^ 

LO = ± — — — I- 



+ O if) . 



(27) 



V3 r(i)r(i) 

This corresponds to a collective excitation of the system - the holographic zero sound mode. 
Its speed is equal to that of hydrodynamic sound, and it has an imaginary part oc g^. In 
the hydrodynamic limit g <^ T the dominant pole is a purely imaginary pole with the 
Fickian diffusion dispersion relation 



iDf + 0{f) 



(28) 



where the diffusion constant is given by |27l. |35| 

d~ 

D{d) = -^Vl + d^2Fi 



3 14 



2' 3' 3' 



:-d' 



(29) 



Between these two extreme temperature limits, the poles of the Green's functions are not 
known analytically. 



20 



35|. The 



These results were generalized to the case of a massive hypermultiplet in 
zero sound mode (1271) persists when the hypermultiplet has a finite mass m, although the 
dispersion relation is altered to 

1/2 



1 - 



I 1 - mV3 



(30) 



r(i)r(i) (i-mV3)^ 

where m = ■m/\/2fl = Mg//i. The real part of the dispersion relation (l30l) was obtained in 
[20I and the attenuation is derived in Appendix |Bl In the hydrodynamic limit, there is again 



a diffusion pole fl28|) whose diffusion constant can be derived via an Einstein relation 



35] 



d-- 

D{d) = —\/d^ + cos^e{l) 



du 



G{u 



cos'^ 6{u)H{u 
sin 6{u) 



l + d\ ^^VV [cose{u)9'{u)] 



G{u) cos 9{u)dd 



+ 



cos^ 6{u) 

where G{u) and H{u) are defined in Appendix Ri 



3 + 



G{u) 



d 



dd 



-sm9{u) 



(31) 
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At finite hyp ermult iplet mass, the D3/D7 system has non-trivial bound states analogous 



to mesons 



36|. These bound states are visible as peaks of the spectral function. When 



39 



46 



48- 



50j. When 



d = 0, such modes exist only for large enough values of m oc Mg/T 
d 7^ 0, it seems that necessary conditions for their existence are large values of m and small 

L-SsJ. This is outside of the regime of our current 



enough values of g/T and ds /T |31 



32 



51 



interest. A full numerical analysis of the quasinormal modes of the theory when T = is 



given m 



54|. 



III. THE HIGH DENSITY REGIMES OF D3/D7 FUNDAMENTAL MATTER 

In this section, we explore the behavior of the collective modes of the D3/D7 system 
in the low temperature regime d ^ 1 {T <^ Anticipating behavior similar to that 
observed in a Landau Fermi-liquid, we expect the system to exhibit the three regimes shown 
in Fig. m and indeed we do find them. We investigate this by computing numerically the 



spectral functions and the poles of the Green's functions using the approach of [3l|, |45 |. 
This approach provides a numerical consistency check (invariance under radial translations) 
which we used to ensure our results were accurate. We always plot the normalized spectral 
functions such as Xzz (t^, q) /2NfNcT'^uj'^: dividing by o)^ has the advantage of reducing the 
high frequency asymptotics ( |26l) to a constant. 



A. The coUisionless quantum regime 

The coUisionless quantum regime corresponds to u and q being in the interval 

< g < 1 , (32) 

see Eq. f|T2l) . Fig. [T]and Table HI We start by fixing a large value of c?, for example d = 10^, 
and computing the longitudinal spectral function at a fixed momentum q in the interval 
fl32|) . e.g. q = 0.4. An isolated peak corresponding to the holographic zero sound mode in 
this low temperature regime is clearly visible in Fig. HI Thus the zero sound mode persists 
at low but non-zero temperatures - contrary to the assertion in that an infinitesimal 
temperature in the field theory will lead to diffusive transport. 

By varying q (while keeping d^ fixed), we can determine the dispersion relation of this 
collective mode at low temperature. Figure |5] shows the real and imaginary parts of the 
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FIG. 4. The holographic zero sound peak in the colhsionless quantum regime. The longitudinal 
spectral function is shown at m = 0, d = 10^, q = 0.4. 



zero sound dispersion relation u = u}{q) at rh = 0, d = 10^. Our numerical results are 

shown as dots and the solid lines correspond to the zero-temperature dispersion relation 

f l27|l . We observe that in the colhsionless quantum regime f l32|) the real part of the zero 

sound dispersion relation shows no noticeable deviation from the zero-temperature result 

( l27l) . The imaginary part shows a close agreement with the zero-temperature result over 

the parameter range ( |32l) . Visible deviations from the zero-temperature result are apparent 

at very small (g < 0.02 ~ d-^/^) and very large (g > 0.7) momenta close to the boundaries 

of the interval f p2|) . Note that although the imaginary part appears to be tending to a 

constant at small q, this is only true up until the crossover to the hydrodynamic regime 

occurs, after which we obtain a diffusive mode. This crossover will be discussed in detail 

later in this section. At large momentum, the system is entering the low temperature, low 

density regime d"^^^ <^ 1 <^ q (i.e. u,q ^ d^ ^ T), where the zero sound mode becomes 

very short-lived as the corresponding pole recedes deep into the complex plane. As shown 

in Fig. ini the zero sound peak in the spectral function gradually disappears in this regime, 

and the spectral function approaches the d = 0, T = result ( !26l) . 

At finite hypermultiplet mass, the results are qualitatively similarjfl The longitudinal 

^ Technically, the massive case is more complicated as it involves solving a pair of coupled differential 

equations, with coefficients that depend upon the numerically-computed embedding function. Due to 

numerical instability, we were not able to obtain accurate results for d > 10^ and masses ouside the 

interval 0.002 < m < 1.68 (at d = 10^). We believe this limitation demonstrates the lack of our numerical 

skills rather than an effect of any physical significance. 
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spectral functions show a lone pealo (similar to the one shown in Fig. H]) for values of uj 
in the range fl32|) . As shown in Fig. [TJ the real part of the dispersion relation is essentially 
identical to the zero-temperature result ( 130|) . and the imaginary part deviates from the 
zero-temperature dependence given in ( 130|) only at the boundaries of the interval ( 132|) . 

It is interesting to note that the zero sound mode exists for all (numerically accessible) 
values of m, including those for which m oc Mq/T ^ 1. Thus the high density and low 
temperature interval in which the T = zero sound mode persists is the same as in the 
massless case (given by the inequality ( |32|) ). i.e. we do not have to take the limit Mq/T — > oo 
which would be the genuine zero temperature limit of the entire system. 



^ Peaks corresponding to "meson" bound states may exist for higher values of w. We have not investigated 
this issue. 
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FIG. 6. The longitudinal spectral function at m = 0, d = 10 . Moving from left to right corresponds 
to increasing momentum: q = 0.4 (black), 0.6 (blue), 1.0 (red), 1.5 (green), 2.0 (orange), 3.0 
(purple). 



B. The collisionless thermal regime 

Having established the existence of a holographic zero sound mode at finite temperature, 
we now investigate the dependence of its velocity and attenuation on temperature, density, 
momentum and mass. The collisionless thermal regime corresponds to excitations with uj 
and q in the interval 

J-^/^ < g < d-^^'^ . (33) 

As discussed in Section IIII A| the imaginary part of the dispersion relation u = u}{q) for 
the zero sound shows significant deviations from the zero temperature ~ behavior in the 
region q < d^^^^. Investigating this further, one can see that in the interval d"^/^ ^ Q ^ d~^^^ 
the imaginary part Imu is essentially independent of q. This is the characteristic behavior 
of the sound attenuation coefficient of a Landau Fermi-liquid in the collisionless thermal 
regime (see Table HTjl . To determine the temperature dependence of the attenuation in this 
regime, in Fig. |8]we plot the logarithm of the ratio lmuj{T)/lmu{0) versus the logarithm 
of d~^^^ ~ T//i at fixed q = 0.01. In the collisionless quantum regime, the attenuation is 
essentially temperature-independent, whereas in the collisionless thermal regime it scales as 
~ T^. The transition occurs at g ~ d~^^^. This is fully compatible with LFL behavior (see 
Table HT] and Fig. |3]). Fig. [8] should be compared to the transition between regions III and 
II in Fig. [31 

At non-zero hypermultiplet mass the results are qualitatively the same and the plots for 
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Log d ~ Log T/d'/2~ Log T/fi 

FIG. 8. The imaginary part of the holographic zero sound dispersion relation at low temperatures 
showing the transition between the collisionless quantum (left) and collisionless thermal (right) 
regimes. The points are the numerical data, the dashed line denotes q = d~^^^ ^ T/fx and the 
solid line is the best-fit straight line for the T/ji ^ q = 0.01 points showing the ~ scaling of 
the attenuation in the collisionless thermal regime. The region corresponding to the hydrodynamic 
regime at even higher temperature T//i S> ^/q = 0.1 is not shown in the figure. 
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FIG. 9. The frequency dependence of the D3/D7 acoustic mode when ?fi = 0, d = 10^ in the 
colhsionless regimes. The dots and crosses are our numerical results, the dashed line is cu = 
irT/d^/'^ ~ T/fi, and the solid line shows the best-fit straight line with gradient a ~ —0.95. The 
points on the left of the left hand plot correspond to the rightmost points on the right hand plot. 



m = 0.76, 1.68 are similar to Fig. [HI We do not show these plots for conciseness. 

Another comparison with Landau Fermi-liquid theory can be made by plotting the atten- 
uation as a function of (real) frequency. As shown in Fig. [2] and Table [Til the frequency 
dependence of the acoustic attenuation, the transition from the colhsionless thermal regime 
to the collisionless quantum regime is characterized by a change of scaling from arg(g) oc 1/u 
to arg(g) oc w at w ~ T//i. We can study this region in the D3/D7 system by fixing the tem- 
perature (at e.g. d = 10^) and tracing the Green's function's pole in the complex momentum 
plane while varying the real frequency u. The results are shown in Fig. [9] for the massless 
case. In the collisionless quantum regime, the points follow a straight line with gradient 
1. In the collisionless thermal regime, for a region where there is a power law dependence 
arg(g) oc a)", the best fit value of a turns out to be a ~ —0.95, which is sufficiently close to 
the LFL value of —1. Fig. [9] can be compared to regions II and III in Fig. [2 As we approach 
the hydrodynamic regime, the power law dependence is lost as expected. We are not able to 
explore the acoustic mode in region I due to limitations of the probe brane approximation, 
as explained in the introduction. The best-fit gradients are the same for non-zero masses 
rh = 0.76, 1.68 and so they are not shown here for brevity. 
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C. The coUisionless-hydrodynamic crossover 



As the temperature is increased further, the zero sound mode becomes less stable, and 
the system enters the hydrodynamic regime characterized by 

< w, g < . (34) 

The zero sound peak in the spectral function broadens and moves to the origin (see Fig. [TOl) . 
The real and imaginary parts of the collective mode's dispersion relation as functions of 
temperature are shown in Fig. [TT] for massless and massive hypermultiplets. The real part 
decreases with increasing temperature until it becomes exactly zero at T = T^ross, and the 
mode ceases to propagate. The magnitude of the imaginary part increases until T = Tcross, 
then decreases again. For T > Tcross, the mode is purely diffusive, approaching at high 
temperatures the known analytic result f l28|) for the hydrodynamic charge density diffusion 
mode. The transition occurs at g ~ rf"^/^, as we will shortly show. 

The dynamics of this coUisionless-hydrodynamic transition are best exhibited by the mo- 
tion of the relevant poles in the complex frequency plane (see Fig. [T2l) . As the temperature is 
raised, the two zero sound poles (corresponding to Rew = if^g) move deeper into the com- 
plex plane, approximately along the circle |t<;| = q/Vs, until they collide on the imaginary 
axis and form two purely imaginary poles. (For the parameters used in Fig. [121 this happens 
i.e. in the region q d ^/^.) One of these new poles recedes quickly 
into the complex plane as the temperature is raised further, while the other approaches the 
real axis and becomes the hydrodynamic diffusion mode ( l28l) at high temperatures. This ex- 
plains the temperature dependence of the collective mode dispersion relation in Fig. [H] and 
the behavior of the density- density spectral function in Fig. [TOl At non-zero hypermultiplet 
mass, the transition is qualitatively similar, although the transition temperature decreases 
slightly with increasing mass {dcross ~ 520, 530 for m = 0.76, 1.68, respectively). 

The phenomenon in which two poles on the imaginary axis collide and generate two 
propagating modes persists at lower densities as well (in fact, it was first observed in Fig. 3(c) 
of [3l| for d = 2 which is the highest value of d studied in that paper), although in the low 
density regime the pole collision occurs deep in the complex plane and the propagating mode 
is short-lived. The pole merger has also been recently observed and correctly identified as 
the hydrodynamic - coUisionless transition involving the zero sound in a 2+1- dimensional 
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FIG. 10. The longitudinal spectral functions for m = 0, q = 0.01. Moving from the top left corner 
to the bottom right corner: each subsequent figure zooms into the area of detail in the previous 
one (note the different scales on the vertical axes). In each figure, moving from the tallest peak to 
the smallest peak corresponds to raising the temperature: d = 10^ (black), 10^ (blue), 10^ (red), 
500 (green), 100 (orange), 10 (purple). 



theory in 23 1. 

In Landau Fermi-liquid theory, the collisionless-hydrodynamic transition occurs when 
~ 1 and g/mfp ~ 1, where r and l^ip are the mean free time and mean free path, 
respectively, of the quasiparticles in the vicinity of the Fermi surface supporting the collective 
mode. Defining the collisionless-hydrodynamic transition in the D3/D7 holographic model as 
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FIG. 11. The temperature dependence of the real and imaginary parts of the dominant cohective 
mode's dispersion relation at g = 0.01 and fh = 0;0.76; 1.68. Dots show the numerical results at 
low T, the crosses show the T = zero sound results (|27p and (j30p . and the solid lines show the 
analytic diffusion result (|28p . 



the event in which the two poles in Fig. [T2]merge, we can cast the corresponding parameters 
in the familiar language of the Landau theory by introducing an effective r = l/lwcrossl and 
^mfp = l/O'cross, and computing their temperature dependence. Assuming a simple power- 
dependence of the form 

Up ~ r oc rf-^/^ (^^) , (35) 
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FIG. 12. The positions of the dominant poles of the density-density correlator in the com- 
plex frequency plane at fh = 0, q = 0.01 as the temperature is changed. The low tem- 
perature limit corresponds to the points furthest from the imaginary axis. As the tempera- 
ture is increased, the points move inwards towards the imaginary axis where they collide to 
form two poles that move up and down the imaginary axis respectively. An animated ver- 
sion of this figure, and the corresponding animations for fh = 0.76, 1.68, are available at 



http : //www . physics . ox . ac. uk/users/Davison/D3D7animations . html 



we expect a plot of log(gcross) or log (| Across I) versus \og{dlloss) to yield a straight line of 
gradient a. These plots are shown in Fig. [13] for the massless case. The numerical results 
are clearly consistent with a simple power law dependence. The slope of the best-fit straight 
line is -2.0 in each case, which leads to the result 

/mfp~r~rfi/3T-2~/iT-^ (36) 



as in a Landau Fermi-liquid 16| - recall table [H For non-zero masses m = 0.76,1.68, the 
best-fit slopes are unchanged and so we do not show the results for brevity. 



IV. SUMMARY AND DISCUSSION 



In this paper, we have investigated the properties of the holographic zero sound mode 
and the fundamental matter diffusion mode in the D3/D7 system at finite temperature and 
high density. Similarly to the case of an ordinary Landau Fermi-liquid, three regimes cor- 
responding to different powers of the parameter T//i <^ 1 can be identified, with the modes 
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FIG. 13. The temperature dependence of the colhsionless-hydrodynamic crossover value of fre- 
quency and momentum for m = 0. The points are our numerical results and the solid lines are the 
best-fit straight lines which both have gradient -2.0. 



having a distinctively different temperature and momentum dependence in each regime. 
In the coUisionless quantum regime, the zero sound attenuation is essentially temperature- 
independent and proportional to the square of the momentum whereas in the coUisionless 
thermal regime, the attenuation is momentum-independent and scales as the square of the 
temperature (see Table HIl Figs. [21 |3] and the corresponding D3/D7 results in Figs. El |9|). 
The crossover transition between the two regimes is clearly visible in Figs. [H |9l In the 
hydrodynamic regime of a Landau Fermi-liquid, the acoustic attenuation is proportional to 
the square of the momentum and scales as ~ Restricted by the probe brane approxi- 

mation, we are unable to study the acoustic colhsionless-hydrodynamic crossover transition 
in the D3/D7 theory directlylJ However, we do observe this transition in the density-density 
correlator as the motion of the lowest-lying poles in the complex frequency plane leads from 
the zero sound-dominated regime to the diffusion-dominated one as the parameter T//i is 
increased (Fig. [T2|) . The transition occurs at g ~ T^/ fi, as if it were supported by Landau 
Fermi-liquid quasiparticles with a lifetime ~ These results do not qualitatively change 

as the mass of the fundamental matter in the field theory is varied, at least in the range 
< m < 1.68 which we have studied numerically. 

The main conclusion of our study is that the holographic zero sound of the D3/D7 
system at finite temperature behaves exactly as the usual Landau zero sound. This is 



See 



55[ for some recent work on going beyond the probe approximation. 
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rather difficult to reconcile, at least within the standard Landau Fermi-liquid paradigm, 
with two other properties of the system: the atypical temperature dependence of the specific 
heat (cy oc /d instead of cy oc TcP^^) and the apparent absence of a singularity when 
|g| = 2qp at zero frequency in the density-density correlator at zero temperature. One 
may add to these apparent discrepancies a finite zero-temperature limit of the entropy 
density 11|. Resolution of the above mentioned puzzles may involve modifications of the 
gravitational background considered, possibly in the spirit of [56]. On the other hand, 
the absence of the singularity in the zero frequency density-density correlator aX 2kp may 
have a simpler explanatioio: assuming the validity of the Luttinger theorem, Ug ^ kp, and 
therefore q 2kp n]^^ ~ since d ^ rin /y/X. This means that the corresponding 

q is large, q ~ A^^^, and is not visible in the correlators in the probe brane limit. It is also 
possible that the ~ behavior of the specific heat is a genuine property of the considered 
microscopic model in the approximation described by the probe brane limit, there is no 
Fermi surface, and we are dealing with a new type of quantum liquid. This point of view, 
initially advocated in is not in contradiction with the findings presented in our paper. 

There are other examples of theories with a dual holographic description which support 
a sound mode at T = 0. It would be interesting to determine whether the results presented 
(q ualitatively) in those cases also. The ffist class of such theories are probe brane 



here hole 
theories 



2ll425l|. Of particular interest are the D4/D8/D8 theory, whose sound mode has 
a different dispersion relation from Landau Fermi-liquid theory at T = 0, and the D3/D3 
theory which is (l+l)-dimensional and hence to which Landau Fermi- liquid theory does not 
apply. A qualitatively different holographic theory is the AdS^ Einstein-Maxwell theory at 
finite density 57-61|- In this theory, non-trivial density- dependent physics is possible 
via the strong coupling of the charge density to the energy density of the field theory. This is 
in contrast to probe theories, where the DBI coupling between the charge density gives rise 
to interesting density-dependent physics, despite the fact that the coupling of this sector 
to the overall energy density is suppressed via the probe limit. These are very different 
mechanisms and hence it would be very interesting to see if the low temperature sound 
mode in the Einstein-Maxwell theory behaves similarly to the corresponding sound mode in 
the D3/D7 theory. 

Finally, it would be very interesting to generalize the recent approach to zero sound 
^ We are indebted to Andrei Parnachev and Pavel Kovtun for valuable discussions on this issue. 
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proposed in [471] to non-zero temperatures. 
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Appendix A: Fluctuations at non-zero hypermultiplet mass 



When the hypermultiplet has a non-zero mass, the bulk fluctuations of the embedding 
scalar and the longitudinal gauge field components are coupled and their solutions are no 
longer independent. In the field theory, this corresponds to mixing of the dual operators. 
An appropriate systematic formalism which we follow was introduced in It will be 

convenient to work with the gauge-invariant variables 



ip [u, uj) = rH(t> {u, u) , E {u, u) = —i [uaziu, u) + qat{u, u)] . 



(Al) 



Note that the definition of E differs from the massless case fl?I]) by a factor of —i. The 
coupled system of equations of motion for the variables E and ip is 



d ( f{u)cos^6 



+ 3iq tan 9A{u)G{u)lp{u, oj 
3 cos^ 9uA{u)9'{u) . 



H{u)E'{u, io) + 4:iqf{u)u'e'{u)A{u)ip'{u, co) 



(A2) 



(is cos^ dH{u) - 



+ 



iqip{u, u) = 0, 



E{u, u) 
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and 



d ( f{u)cos^6 



du yu^/G{u)D{u,u}) 
- 4iqf{u)u^e'{u)A{u)E'{u,0) - 3ta.nee'{u)G{u) {u"^ - f f{u)) ip{u,oo) 
3 cos^e sin ef{u)A{u)^/G{'t 



(A3) 



+ 



D{u, u) 
cos^ euA{u)e'{u)d 



iqE'{u, u) 



iqh[u, CO) H (p [u, 00) 



u 



uDiu, u) 



+ 



cos^ edz ju'^ -ff{u)B{u)) 



gcos^sin^ eA{uf^/G{v) 



D{u, oo) 



3(3 cos3^-2cos^) JG{u) , 
+ = 0' 



where the coefficients are given by 

A[u) = ^^, B{u) = \-^u^A{uf , G{u) = \^^u^f{u)e'{uf -^u''A{uf 



H{u) = 1 + 4f{u)u^9'{uY , D{u, CO) = u!^H{u) - ff{u)G{u), 

and primes denote derivatives with respect to u. 

The relevant part of the off-shell action quadratic in longitudinal fluctuations is 



c(2) 
"-"long. 



doodq I cos^ 9 f{u)H{u) -, 



(27r)M ^/G{^D{u,i:o) 



E'{u,-u)E'{u,co) 



(A4) 



1C0s3^/(m)W(m)A(m). , ,^ , 

^ ^ ^ ^ ^ -iqE'{u,-uo)ip'{u,u) 



+ 



^/G{v)D{u,Co) 
cos'efju) {co'-q'f{u)B{u)) 

u^^yG{u)D{u, CO) 
6 cos^ 6 sin ef{u)A{u) ^/G(^] 

6 cos2 e sin ef{u)e'{u) ^/G{^ (a;2 - ff{u)) 



uD{u, uj) 
+ non-derivative terms 



iq ip{u, —uo) E'{u, u) 

{p{u, —u) (p'{u, u) 



where N = N fTo'jVss. The equations of motion and the action are written in the form that 
allow one to apply the recipes of 3l|] directly. To obtain the retarded Green's functions, we 



solve the coupled system of equations (lA2p - (1A3P with incoming wave boundary conditions 
at the horizon and combine with the appropriate factors from the action (]A4|) as described 



m 



3l|. 
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Appendix B: Holographic zero sound attenuation at finite hypermultiplet mass and 
zero temperature 



In this appendix we derive the formula for the zero sound attenuation at finite hyper- 
multiplet mass following the approach of Ref . {20!] . Note that our notation is different from 
that used there. At zero temperature and non-zero mass, it will be convenient to use the 
coordinate system in which the background metric takes the form 



— {-dt^ + dx^) + — [dp^ + p^ds% + ci7^2 + Tl^dct)'^) 



where = + TZ^. The background DBI action is given by 

/■oo 

Sfund. = -N dpd^x p' ^Ji + n'ipY-A^ipf, 

Jo 

where 71 is the embedding coordinate of the D7-branes. The background gauge field and 
embedding coordinate are determined by two conserved charges c and d via 

i?2 



Aip) = —d{d'-c^) -^B 



n{p) 



6 



c{d^-^y^'B 
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1" 


[p6 


+ d2i?12 _ c2i?12 
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3 




P' 
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1" 


[p6 


+ c/2i?12 _ c2i?12 ' 


6' 


3 



where d is the density of fundamental matter and B is the incomplete beta function. The 
mass and chemical potential are given by the asymptotic values 



At (p ^ 00) = vhP, 



and are related to the constants d and c via 



V2 



d = rlrP {~p' - rhy2) 



where 



7 



1 1 
6' 3 



-3 



The equations of motion for the fluctuations of the embedding scalar 6TZ and the fluctuations 
of the gauge-invariant combination 6 A = R^uaz + R^qat are couplecj^. In the near-horizon 
limit p — 0, the solutions take the form 

5A = Ap e^*^/", 5n = Bp e^'^^', (Bl) 

— n 

^ They are given by equations (A. 11) of [20'|, with the replacements E — > SA, ^ — > STZ, r ^ p, d ^ dR^, 
cR^, R, R^n. 
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where f2 = B?l>j\/\ — c^/d? and the upper sign in the exponent corresponds to the ingoing 
boundary condition. In the small frequency limit (fi/p <^ 1), the solutions (IBip can be 
expanded as 

5 A = ±inA + Ap+---, 6n = ±iVLB + Bp + -- - . (B2) 

Alternatively, taking first the small frequency limit of the equations of motion, we obtain 
the following solutions near the boundary 

5A=Co + (l/p2) , 5n = Co + (l/p2) , 

whereas near the horizon the corresponding solutions are 



5A = Co + hCi + 
6n = Co + hCi + 

where the coefficients are given by 



aiCip + a2C2P, 
aiCip + a2C2P, 



(B3) 



ai 



hi 



hi 



r (i) r (I) [(3c2 - <e) + 3 (d^ - c^) u'] 



02 



cdq^ 



i?2 (d2_c2)2 

2cdq-T (I) r (I) 



r (1) {d^ 



cd 



r (I) id' 

rf2 



^cdT (i) r (I) 

4 ; 

-2^3 



02 



3 ' 
^2^2 



i?6 (rf2 _ c2 

r (I) r (I) (3^2 _ ,2) 



/24r (I) (rf2 



4 

-^2^3 



i?4r (i) (c/2 _ 

The poles of the retarded Green's function correspond to the quasinormal modes of the 
system - these are the solutions which obey ingoing boundary conditions near the horizon 
and vanish near the boundary. Such solutions exist when it is possible to match the expansion 
(lB2p (with the upper sign) onto the expansion (1B3P with Co = Cq = 0. Performing this 
matching order-by-order in p, we find that the coefficients Ci and C2 obey the equation 




iVlai — bi iVLa2 — &2 
iVtdi — hi iVLd2 — &2 

A non-trivial solution exists provided that 



iVtai — hi iQa2 — 62 
iQdi — hi iQd2 — 62 



(B4) 



(B5) 
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By substituting an expansion of the form 1^ = 00 + Ciq + C2q^ + . . ., we find that the above 
determinant vanishes when u is given by the dispersion relation fl30l) . Note that this method 
only finds the lowest energy quasinormal modes, as it involves taking the small uj limit. 
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